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Abstract 

We consider recovering analytically the (generally complex) parameters /x, Mi and 
M2 of the gaugino and Higgsino Lagrangian, from appropriate physical input in the 
chargino and neutralino sectors. For given tan/5, we obtain very simple analytic 
solutions for M2, |//|, ^r(7[/i] in the chargino sector and a twofold \Mi\, Arg[Mi\ 
analytic solution in the neutralino sector, assuming two chargino, two neutralino 
masses, and one of the chargino mixing angles as physical input. The twofold am- 
biguity in the neutralino parameters reconstruction may be essentially resolved by 
measuring the e"'"e~ X1X2 pi'oduction cross-section at future linear collider ener- 
gies, which we study explicitly with the phase dependences. Some salient features 
and specific properties of this complex case gaugino "spectrum inversion" are illus- 
trated and compared with the similar inversion in the real case. In particular, our 
algorithms exhibit in a direct and transparent way the non-trivial theoretical cor- 
relation among the chargino and neutralino parameters, and the resulting allowed 
domains when only a subset of the required physical input masses and production 
cross-sections is known. 
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1 Introduction 



In the Minimal Supersymmetric extension of the Standard Model (MSSM) [|I], 0, with- 
out additional theoretical assumptions, the soft-supersymmetry breaking part of the 
Lagrangian|^ involves a large number of arbitrary parameters. Moreover, many of these 
parameters may be in general complex, adding new sources of CP violation with respect 
to the standard model 0, Apart from the many possible phases in the flavour sec- 
tor which are mostly severely constrained |P, ^ |^, the other "flavour-blind" possible 
phases, appearing in particular in the gaugino and/or Higgs sector, are essentially con- 
strained P] by the electric dipole moment of the electron P] and the neutron [l^. As 



recently emphasized p], |T2|, large phases in this sector are not excluded at present, al- 



though it may be considered as unnatural because it generally requires specific cancella- 
tions among the different contributing phases [1T3|, p^ . Irrespective of a realistic theoretical 



scenario implying such phases, the latter lead, from a more pragmatic point of view, to 
potentially drastic changes in the phenomenology of the Higgs and gaugino sectors of 



the MSSM, therefore affecting the already challenging reconstruction from data [[15|-||19 
of the structure of the SUSY and soft-SUSY breaking Lagrangian. Typically, assuming 
that a certain amount of the SUSY partners may be discovered and some of their masses, 
production cross-sections and decay widths measured, additional questions arise in the 
presence of non-zero phases. Even without direct measurement of manifestly CP-violating 
observables, one may for example try to determine whether a set of available data neces- 
sarily implies the existence of non-zero phases, or conversely whether the data is consistent 
with the underlying MSSM basic Lagrangian with phases. 



In a previous paper [|T^, we have given a specific strategy and algorithm to de- 



diagonalize or "invert" the parameters in the ino sector, i.e. obtaining some of the 
MSSM Lagrangian parameters in direct analytic form in terms of physical parameters 



(see also ref. |]2^ for an early such attempt). However, our construction was restricted 
to the case of real parameters. In the present paper, we generalize our procedure to 
the complex parameter case. While the techniques used are quite similar, especially for 
the less straightforward 4x4 neutralino mass matrix inversion, the choice of the most 
appropriate input is obviously different, since the number of input/output parameters is 
necessarily larger when phases are present. In addition, some peculiarities of the complex 
case require a slightly different approach. As we shall see, it turns out that the complex 
case inversion is much more constrained than one may think at first: strong correlations 
between the chargino and neutralino physical masses occur solely from the structure of 
the MSSM Lagrangian, despite the additional phase freedom, and we will illustrate how 
this somewhat unconventional spectrum inversion algorithm can exhibit in a more direct 
and clearer way such non-trivial correlations. For instance, it becomes straightforward 
in that way to exhibit allowed domains for the neutralino masses, once some of the pa- 
rameters of the chargino sector are determined (or vice- versa), a task that would usually 
require rather tedious systematic scanning over the basic Lagrangian parameters with 
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painstaking comparisons of physical output results. 

For definiteness we shall only consider here the reconstruction of the chargino and 
neutralino sector parameters with phases, letting aside e.g. the possible implications of 
non trivial phases on the MSSM Higgs sector phenomenology, which is somehow decou- 



pled from our analysis and has been investigated recently in refs. [21, 22]. Also, as in 



ref. []T^, we emphasize that we are not aiming at a new specific MSSM parameters re- 
construction method from experimental observables. Clearly an algebraic approach alone 
cannot replace more complete studies of the reconstruction of the MSSM parameters from 
data. Our algorithm should rather be viewed as an efficient tool, to complement any more 
systematic analysis of the data. 

Our brute algorithm gives analytical solutions in the unconstrained MSSM (i.e. no uni- 
versality of the gaugino mass terms is assumed), for /i = |/i|e**'', M2, and Mi = |Mi|e***^i , 
using as input two chargino masses M^^, M^^, two (arbitrary) neutralino masses M^^, 
M]s[. , and one chargino mixing angle, plus tan (5 As will be explained, we can choose 
a phase convention such that M2 is real, so that the counting of physical input versus 
output parameters is consistent. Obviously, any of the latter input, if not known, may be 
scanned over some specified range to determine e.g. the resulting allowed domain in the 
output parameters /x. Mi, M2. As a by-product, we obtain the neutralino diagonalization 
matrix elements -and therefore the neutralino physical couplings - as direct functions of 
the masses, and we have recalculated the e^e~ — > x^X^j production cross-section in this 
framework with complex /i. Mi. The neutralino 1,2 pair production cross-section has 
been widely studied as a very promising detection process at LEP2 ||2^, |2^ or at a 



future e"'"e~ linear collider |2^, The next- to-lightest M^^ mass could be measured 
with good accuracy, and the lightest (LSP) M^^ as well from chargino decay product [p^ . 
As we shall illustrate, the neutralino pair production measurement should also play an 
important role in resolving an intrinsic ambiguity in the reconstruction of the Mi param- 
eter from our algorithm. One should note, however, that prior to our analysis, the precise 
reconstruction of the chargino masses and mixing angles from actual data will be non 



trivial. This issue was recently investigated in details in ref. |T8|, 0, with another way of 
reconstructing the /i and M2 parameters by considering the chargino pairs which can be 
produced at a future e^e~ linear collider and studying theirs spin correlation properties. 
Accordingly, our analysis is very complementary to theirs, in particular for the neutralino 
sector. 

In section 2, we briefiy recall the general parameterization of the gaugino and Higgsino 
sector, and the physically relevant phases fixing our conventions. In section 3, we discuss 
the complex chargino and neutralino de-diagonalization algorithm and discuss uniqueness 
conditions. Typical illustrations of the reconstruction of the parameters are given, and 



in particular the kind of non-trivial correlations that one can glean in the case of partial 



knowledge of the input. In section 4 we give analytical expressions for the e^e — > Xi'x? 



^This algorithm is available as a fortran code upon request to the authors. 
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production cross-section in this general complex /x, Mi case and illustrate how it can re- 
solve the two- fold ambiguity for Mi, when determined from the neutralino masses alone. 
Finally section 5 gives some conclusions and outlook. 

2 Phases in the MSSM gaugino Lagrangian 

The relevant Gaugino/Higgsino part of the unconstrained MSSM Lagrangian can be found 
in many different places, see e.g. refs[|l], |^. More specifically the soft SUSY-breaking Ma- 
jorana mass terms: 

C„no = -^BB - ^WW, - ^G'^Ga + h.c. (2.1) 

have to be supplemented by the supersymmetric term originating from the quadratic 
part of the superpotential and contributing to the Higgsino terms. Writing explicitly 
for illustration only the terms contributing to neutralinos (after electroweak symmetry 
breaking) : 

C-neutralino = mzC^j siu PW^Hu - mzCyj COS PW^Hd 

+mzSw COS (3BHd - mzs^ sin /3BHu + ^HuHd + h.c. 

(2.2) 

with Sw = sin 9w, Cw = cos 9w, and tan/3 = Vu/vd the ratio of the two Higgs vacuum 
expectation values, Vu^d = {Hu,d), it is now easy to see from Eqs. (|2.1| , pl^) that only the 
relative phases of say (M2,Mi) and {M2, fi) are physically relevant. Indeed any phase 
change of M2 in Eq. (|2.1| ) can always be absorbed by a phase change of the W field. The 
latter, however, fixes uniquely the phase change of Hd, Hu and B in Eq.( p.2D since we 
assume all vev's and gauge couplings to be real, in such a way that the phases of the 
combinations M2/M1 and /iM| remain unchanged. We can thus choose, without loss of 
generality, M2 > and real, and Mi and to have arbitrary phases, which fixes our 
phase convention. Even though no symmetry arguments were invoked here for Cneutraiino, 
one should keep in mind that Peccei-Quinn and R-symmetries are implicitly needed to 
assure that phase rotations of the fields, necessary for a given phase convention, preserve 
automatically the reality of quantities appearing in the PQ- and R-symmetry preserving 
sector of the MSSM, such as quark and lepton masses and gauge couplings.0 

3 Extracting the parameters /i, M2 and Mi 



^Note that the assumption of no phases for the Higgs vev's is the only convention here that affects 
the phases of the remaining soft susy sectors of the MSSM. For a clear discussion on the relevant MSSM 
phases from C/(l) symmetry arguments see, e.g., refs. 
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3.1 Chargino sector 

In the chargino sector, and with the phase convention discussed above, the mixing mass 
matrix reads from Eqs. ( |2.1| , |2.2| ) 

The squared chargino masses are obtained as the eigenvalues of the MqMc matrix: 



^xi,2 = ^[^2 + l/^r + 2m2^ T A ] (3.2) 



where (A = M^^ - M^J 



A = ^(M| - |/i|2)2 + 4m^cos2 2/5 + 4m2^(M| + \fi\^) + Sm^^Msl/x] sin 2/5 cos . 

(3.3) 



Note that the diagonahzation of the non-symmetric matrix ( p.l|) is performed via two 
unitary matrices, involving mixing angles 



UlMcU],' = Mc,d^ag (3.4) 



with 



U = ( e *^^sin0L\ 

^ y —e'^^^silKpL COS0L J 

^ ~ 1^ e^^2 j 1^ -e^/5« sin cos0/j^ j 

where the phases (3l,r and 71,2 depend on the parameters M2, (and tan/5). Their 

explicit expressions, that we shall not need in our subsequent analysis, may be found e.g. 
in refs. [0, 01- 

From the relations UIMcMI{UI)~^ = UrMIMcU]^^ = diag{M^^, M^J one obtains 

M|- |/ip-2m^^ cos 2/5 
cos20L = — — (3.6) 



2mwJMl + |/i|2 + (M| - |/i|2)cos2/5 + 2M2|/i| sin 2/5 cos 
sin 2(f)L = 

and obtained by similar expressions in which sin/5 ^ cos/5. 

Accordingly, in the complex case there are two independent parameters among (pi, (j)^ 
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For convenience we adopt notations similar to those of ref. |l9l. Note, however, an overall difference 



of signs of (j)L, 4>R angles in eqs. ( |3.6[) , which should be due to a different convention for the ordering 



of the two chargino mass eigenvalues. Our convention is such that M^Mq = diag{M'^^, M^^) (with 



and tan /?, in addition to the two physical chargino masses. The mixing angles 0^ and (pR 
enter in different combinations in the ZxfxJ ^Xi^e couplings, such that, as worked 
out in details in ref. |19i, precise measurements of the e^e~ — > xtXj cross-section 



and final state polarization properties should allow a non-ambiguous reconstruction of the 
two chargino masses, mixing angles 0i and 0/?, and even tan/3. We refer to their analy- 
sis for further details. For convenience, in our subsequent analysis, we rather choose to 
have tan (3 as input, having in mind that tan (3 maybe also determined from another sector. 



Therefore, assuming as input tan/?, the two chargino masses and one mixing angle 
(say (pL for definiteness) , one can derive straightforwardly from eqs. ( |3.6|) "inverted" 
expressions for $^ and M2 directly in terms of physical parameters. These read 



2^ 



2m 



cos 2/? - (+)Acos(20i)) 



(where S = M^^ + M^J and 



cos($^ 



P = M2 



1/2 



M^^M^2^-(P-m^ sin 2/3)- 
2m^ P sin 2/3 



(A cos 2(pL + 2m^ cos 2(3f - (S - 2m^) 



1/2 



(3.7) 



(3. 



Actually, since only cos($^) 
to a twofold ambiguity: $^ 
manifestly CP violating observables|19 



enters the mass formulas ( p.2[) , $^ is determined by ( |3.8|) up 
^ 27r — To resolve it requires a priori the measurement of 
Note however that the only effect of this ambigu- 
ity on our chargino, neutralino inversion algorithm is accordingly $jv/i ^ —^Ah, without 
affecting \Mi\ or the e'^e~ — ^ x'iX'j production cross-section, as will be illustrated in 
subsequent sections. 



Now, what is not immediately transparent from eqs. ( p.7| , p.8|) is that they already 
exhibit relatively strong constraints among the physical parameters. Namely, arbitrary 
values of M^^, M^^ and < 0^ < 7r/2, are not at all guaranteed to be consistent with the 
obvious constraints M2, > 0, and | cos^^l < 1. More precisely, one arrives, after some 
straightforward algebra, to the two consistency constraints 



^A^ cos2(20i) + 4m^A cos 2/3 cos(20i) + 4m^S - A^ 
< 16m^ sin2(2/3) (A cos(20l) + 2m^(l + cos 2/?) - S) 
X (A cos(20i) - 2m^(l - cos 2/3) + S) < 



^m^ sin2(2/3) 



Acos(20i) + 2m^(l + cos2/3) -S < 



(3.9) 
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In other words, only those values of M^^, M^^ ^.nd cos20i which satisfy eqs. ( |3.9| ) 
are consistent within the MSSM, which for some masses may give rather restricted ranges 
for cos20i. These algebraic constraints mean that, given and M2, one cannot reach 
arbitrary values for the chargino masses despite the phase freedom, the magnitudes of 
these masses being essentially determined by |//| and M2. Obviously, the actual chargino 
pair production will correspond to only one consistent (pi value, but ( ^.91 ) also gives in a 
straightforward way some other useful information, like for instance the allowed (M^j, tpi) 
domain if only the lightest chargino mass is known. Illustrations of such purely theoretical 
consistency constraints will be given in the next sections. 



3.2 Neutralino sector 

Let us now turn to the de-diagonalization of the neutralino sector. The question we want 
to answer analytically here, is how to determine in general \Mi\ and (as well as 
two remaining neutralino masses), for given M2, tan/? and two arbitrary input 

neutralino masses M^r^, MjVjQ- At this stage /i and M2 may be either obtained from the 
previous procedure, eqs. ( |3.7| -|3^), or known by any other mean. Clearly, once the latter 



Ml reconstruction algorithm is obtained, it can easily be combined with the one of the 
previous section to determine Mi,M2 and /z directly from input values of the physical 
parameters M^^, M^j; 4'Li M^-^ and Mjy^, as will be illustrated in section 3.4. 
The neutralino mass matrix with the relevant phases reads 

/ |Mi|e**"i — m zsw cos P mzSwsmP \ 

M2 nizcw cos (3 —mzcw sin (3 

—nizSw cos P rrizCw cos P — |/i|e**'' 

V m^siysin/? —mzCw sin P — |/i|e**'' / 



M 



(3.10) 



Since M is now complex and symmetric but not hermitian, it cannot generally be di- 
agonalized through a similarity transformation (see eq.( [4.1|) ). In order to proceed as in 
ref[|l3, we should rather consider the hermitian matrix M^M, and use the following four 



invariants under a similarity transformation: 
T = TrM^M = ^k'^ 

4 = 1,4 

{TrM^Mf Tr{{M^Mf) 



2 2 



{TvM^Mf TtM^M TT{{M^Mf) TrjiM^Mf) _ , . , 

7y 1 n — 2^ ^^-^ Ni-^^-^ N/^-^ . 



D = DetM^M = U Ml : 

1 = 1,4 ' 

^In what follows we denote the two input neutralino masses by Mjvi , M^f^ for definiteness and assume 
them to be the lightest and next to lightest. It should be clear, however, that any two masses among the 
four will can be equivalently used as input. 



(3.11) 
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where now the M|r. are the (squared) mass eigenvalues satisfying the eigenvalue equation 

det{M^ M - Mil) (3.12) 
with coefficients of (M^)* given by ( p.ll|) for i = 1, ..4. 

After some lengthy but straightforward algebra, we obtain the explicit expressions 



for (|3.11| ) as a function of the various relevant parameters. Following then a procedure 



similar to the one discussed in the appendix of ref . [|T^] , but with the obvious replacements 



M —>■ M'^M and M^. — > M^^, we can solve this over constrained system analytically to 
give I Ml I, and the two remaining physical neutralino masses. It is convenient in a 
first stage to give the form of the equations controlling Re[Mi], Im[Mi\: 

aiRe[Mi] + biIm[Mi] + ci\Mi\'^ = di 



a2Re[Mi] + h2lm[Mi] + CajMij^ = d2 



(3.13) 



where the coefficients Oj, 6j, q and di are functions of the parameters M2, Mat^, 

Mjvj and tan/3. Once i?e[Mi] and Jm[Mi] have been determined from (|3.13|) , the remain- 
ing unknown (squared) neutralino masses M^q.^^, M^r^ are obtained as the two solutions of 
a quadratic equation: 

1 \t - - Ml^ T ({T - Ml^ - Ml? - AD/{M^,M^;^'^'''^ 



(3.14) 

where the trace T and determinant D are given as functions of /x, Mi, M2 from eqs. (13.11]). 



The algebraic expressions for |Mi|, $4/^ are not particularly telling. We thus refrain from 
giving here more than their generic form explicitly (their fortran encoding does not take, 
though, more than a couple of dozens of lines). 

Rather, the more interesting characteristic features of this neutralino inversion is that, 
quite similarly to the simpler chargino case discussed above, non-trivial constraints among 
the physical masses arise, simply because the system ( p.l3| ) cannot always have a solution 
for any arbitrary Mjvi, Mjq^ input. More precisely, the system is (conditionally) either 
over constrained and has no consistent solution; or has a twofold solution. This again 
simply reflects that not all possible M^Vi, Matj values can be consistently reached when 
taking into account that | cos$MiI, I sin^jvfj < 1- (For instance it is clear that the $Mi 
phase freedom cannot compensate for arbitrary values of M2 or |Mi|, which more 
essentially determine the magnitude of the two neutralino masses). In the specific form 
taken by the solution of the system ( |3.13| ) this issue is entirely driven by the positivity 



condition of the radicant A^v, entering the quadratic solution for Re[Mi]: 

Aat > (3.15) 

Re[M,] = ^±y^ (3.16) 
2a 
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Im[Mi] = jRe[Mi] + S (3.17) 

where A^,a,P,j and S have a comphcated dependence on M2, M^Vi, M7V2 

and tan/3, which we do not give here. [The quadratic structure of the solutions is due to 
the special form of (|3.13|) where the non-linear terms appear only in |Mip, otherwise one 
would have expected a quartic structure.] Accordingly, whenever the system is consistent, 
i.e. eq.( 3.15| ) satisfied, there are two solutions. It may seem at first sight that this ambigu- 
ity occurs only as an artifact of "squaring" the mass matrix in our method. This is not so, 
for instance even if one considers both Mi solutions afterwards in the actual (unsquared) 
mass matrix ( p.lO|) , the ambiguity is not resolved but will be simply reflected in some 
phases of the mass eigenvalues, or equivalently in the couplings of the corresponding neu- 
tralino species. It follows that this twofold ambiguity on Mi, intrinsic to the case where 
only two neutralino masses are input, may in fact be resolved by a study of the corre- 
sponding neutralino pair production cross-section which we consider in details in section 4. 



3.3 Reconstruction of Mi, an illustration 

To put more flesh on the previous algebraic construction, we shall illustrate here some 
scenarios with typical input choices, and study the consequences on the relevant param- 
eter reconstruction. Let us first consider for simplicity the pure neutralino parameter 
|Mi|, <l>Mi reconstruction, thus taking \fi\, $^ and M2 as input. Even though this is 
not a reconstruction of Mi directly from physical parameters, this partial algorithm has 
some interest before we consider a complete reconstruction from physical chargino and 
neutralino masses. 

In fig. we chose typical input values |/i| = 100 GeV, M2 = 120 GeV (and tan/5 = 2), 
such that one of the charginos is relatively light, M^^ ^ 50-120 GeV, very roughly within 
the LEP2 reach or exclusion range Q We plot the resulting solutions for Mi {\Mi\ in 
fig la (top) and in fig. lb (bottom)), as a function of the phase and for fixed 
MjVj = 40 GeV and M^^ = 80 GeV input. (Note that in fig. lb, the discontinuity in one 
of the solution is not physical, being only an artifact of the definition of the phase 
modulo 27T. Our convention is such that — tt < < tt). The shape of the different 
plots is rather generic, and several characteristic features can be noticed. First, there is 
the rather important central zone where no solution for Mi exists, corresponding to the 
consistency constraint ( p.l5|) not being fulfilled. For the input choice in fig. ||, this no- 



solution zone covers about one third of the possible $^ values, which is not so restrictive, 
but it should be noted that a relatively moderate change in the value of either M2/\fi\, 
(tan P), or M^^^/M^^ may easily reduce considerably the consistent region, and quite often 
no consistent Mi can be found for any $^ within a relatively important domain of the 

^One should, however, keep in mind that the existing bounds||^, ^ |2^, on top of being model 
dependent, have been derived assuming real parameters. 
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input parameter space. This will be illustrated more quantitatively below with "scan- 
ning" plots. One can already see the sensitivity of the consistency relation ( p.l5| ) to those 
changes in fig 2a, 2b, where a larger tan/3 and a lower Mjvj neutralino mass have been 
chosen. 

In the range where Mi solutions do exist, the plots in fig. la,b and 2a,b clearly 
exhibit the twofold ambiguity on Mi (except obviously at the $^ value (s) such that 
Atv = 0, the boarder of the no solution zone, where Mi is unique and the two distinct 
curves should meet.) For the input parameter choice in figs. 1 and 2, one observes that 
the two solutions are much more distinct for $Mi; than for |Mi|: this is rather generic, 
although larger differences in the two |Mi| solutions could occur, for instance for very 
different M^v^, M^^ values. In fact, as is clear from the figures, the order of magnitude of 
|Mi| is essentially determined by the values of M2 (or M^J ~ O (100 GeV), while the 
$^ variation here has a rather moderate infiuence on |Mi|. For example in fig. la, |Mi| 
varies as a function of $^ roughly between a minimal (~ M^Vi) and a maximal (~ M^vg) 
value, but that property is somewhat accidentally due to the relatively close values of 
and M2 chosen here. When M2/\fi\ is larger (or smaller), we observed that the consistent 
values of Mi are very close to the lightest M^-^ input values, although it is hard to derive 
a simple generic behaviour. 



3.4 Complete gaugino inversion, an illustration 

We will now illustrate the merging of the two separate algorithms described in previous 
sub-sections 3.1 and 3.2 for the chargino and neutralino sectors. The outcome is thus 
the values of M2, |Mi| and $Mi as direct expressions of the chargino masses and 

mixing angle (pi, and of the two arbitrary neutralino masses (plus tan/3). 
In fig. 3a and 3b are plotted the output modulus and phases respectively, as functions of 
the physical chargino mixing angle 0^ for typical input chargino and neutralino masses 
(see figure captions), the lightest being roughly close to present experimental exclusion 
limits (although again, the exclusion analysis performed assuming all parameters real may 
not apply here). The non-trivial correlations among the chargino and neutralino masses 
appear now explicitly. First, and rather generically, the domain where both sectors can 
consistently exist is quite narrow as a function of (pi (the right boarder zones in fig. 
3a,b are first eliminated from inconsistency with the relations ( |3.9| ), while the condition 
Atv > eliminates most of the remaining (p^ domain, leaving the relatively narrow Mi 
solution zones). Again, relatively moderate changes in the input mass values may easily 
result in narrower or even empty solution zones. For instance, changing M-^^ only from 
80 to 100 GeV, for the same values of the other input parameters, gives no consistent 
solutions. Another example is fig. 4 which illustrates the same reconstruction as in fig. 
3 but for a slightly higher tan/? value. In summary, plots such as the ones in fig. 3 
and 4 show that the inversion provides in a direct way the kind of correlations one is 
expecting among the neutralino masses, once we know the chargino spectrum (or vice- 
versa). We could obviously have illustrated the same kind of information by choosing for 
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instance one of the physical masses to vary, instead of 0/,: the plots would be different 
but would qualitatively reflect the same strong correlations. This will be illustrated more 
systematically in the next sub-section. 



3.5 Reconstruction from partial input 

We consider now a less optimistic situation where only a partial knowledge of the physical 
input parameters is assumed, and illustrate the kind of information that can be retrieved 
in this case. Accordingly in figs. 5-8 some input masses are fixed while other and/or 
mixing angles are randomly scanned within a reasonable range for the masses (see figure 
captions). First, in fig. 5 we fix the two chargino masses and tan/5, and vary (pi between 
(pi ^ 0.37 and (pi ^ 0.51 (rad), which corresponds to the first consistent solution zone in 
fig. 3a,b. Condition (|3.15|) implies the pattern of correlation among the physical neutralino 
masses, represented by the dotted regions ("butterflies") in the figure. More precisely, a 
given set of the four neutralino masses is consistent only if any pair (M^v^, Mjq.) of these 
masses corresponds to a point lying on one of the "butterflies". Definite consequences 
follow from this requirement such as the fact that each of the three allowed branches 
(along say the y-axis) can host only one pair (Ma?-, M^.). This allows to make qualitative 
statements about the neutralino spectrum relative to that of the charginos. For instance, 
with the specific input chosen in fig. 5, close inspection shows that there should always be 
one (and only one) neutralino heavier than the heaviest chargino, and that — Mat^ 
gets smaller when M^^ — M^^ gets larger, and vice-versa (M^r^ is the LSP and M^^ the 
heaviest neutralino in our conventions). If tan/3 increases, the dotted butterflies are sim- 
ply moving up or down along the diagonal M^r- = M^r^ line, as illustrated for the choice of 
parameters in fig. 6: thus the (anti) correlation property between the lightest and heaviest 
mass splitting remains, but for larger tan j3 values there is room for the four neutralinos 
being lighter than M^^. Finally, note that for the input parameter choice of fig. 5 and 6, 
there is no consistent solution when tan /3 > 8. 

We stress again that these correlations in fig. 5 and 6 are only due to the theoretical 
consistency relation ( p.l5|) , and are simply obtained from scanning over the two input 



neutralino mass ranges, calculating ( p.l5| ) only once for each M^^., M^., (and eventually 
(pi) input. On top of this, one could easily add any direct experimental constraints on 
the LSP and other gaugino masses |^ Q etc (which we however refraine from investi- 
gating here, since it is difficult at present to infer which experimental mass limits remain 
valid when non-zero phases are assumed). From a practical point of view, it may be use- 
ful to compare how figs. 5 and 6 could be obtained from a more conventional procedure. 
Accordingly, one would have to scan over the five-parameter space M2, \Mi\, $Mi, 

with tan/3 = 2 (tan/5 = 5), diagonalizing the complete chargino and neutralino matrices 
for each input point, and finally selecting only those output corresponding to the required 
physical chargino and neutralino masses. Even if our method is not intended to replace 
conventional procedures of e.g. fitting to the data, it can be in many practical situations 
an appreciable "theoretical consistency" guide with a substantial gain in simulation time. 
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Next, we consider another application of the gaugino spectrum algebraic inversion, 
where we try to illustrate the kind of minimal data needed to e.g. draw definite conclusions 
about non-zero phases. Figs. 7 and 8 give typical "contour" plots in the slice 
of the basic parameter space, when assuming that only the lightest chargino and the two 
neutralino masses are fixed, while the heavy chargino mass is basically unknown, and 
also tan/? varies within some range. Figs. 7 and 8 illustrate two different choices for 
the chargino and the two neutralino masses in this scenario (see figure caption). For the 
input chosen here, one can see that the "data" definitely imply large non-zero and 
this statement does not depend on the specific value of M^^, as shown by comparing 
fig. 7 and 8 input. However, when increasing tan/5 or the splitting between the two 
neutralino masses, more and more points appear with $^ close to zero. Qualitatively, 
one can summarize this behaviour by noting that there exists a well-defined correlated 
region in the chargino-neutralino mass parameters, such that for small enough tan /3, 
is necessarily large (the precise value of tan/? depending on the precise values of the 
neutralino mass splitting). It should be thus interesting to simulate this scenario in a 
more systematic way, taking into account also independent constraints on the phases 
$Afi m^iHI, on tan/3, etc, which is however beyond the scope of the present paper. 



4 The e+e X?Xj cross-section with phases 

In this section, we recalculate and illustrate the neutralino pair production in the complex 
parameter case, in e~^e~ collisions. Our motivation is twofold: 

i) once non-zero Mi, /i phases are assumed in the gaugino sector, all predictions for 
chargino, neutralino production processes will be changed, and it is crucial to analyze the 
relevant production cross-section with the general phase dependence. This was done for 
chargino pair production in reffl^, |19| precisely to extract the chargino sector complex 
parameters. As for the neutralino pair production in e~^e~ collisions, the cross-section 
with non-trivial phases was calculated in the past p9[ and CP-violation effects in the 



neutralino pair production and decay was also considered in some details 

ii) As was stated previously, the reconstruction algorithm, using only two neutralino 
mass input , leads to a twofold ambiguity in the resulting Mi values. We shall illustrate, 
however, that the above-mentioned Mi ambiguity can be in most cases resolved from the 
corresponding neutralino 1,2 production total cross- sect ion, if the latter is measured with 
a reasonable accuracy. As a by-product of the inversion procedure, we can also obtain 
direct correlations between the chargino and neutralino pair production cross-sections. 
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4.1 Analytical expressions 



The neutralino i,j pair production cross-section in e^e~ collision at tree-level proceeds 
through s channel Z, 7 and t(u)-channel sneutrino Ue exchange. The cross-section for 
the real parameter case was calculated in ref . pT|, E^, and for the complex parameter 



case in refs [|^. Here we have recalculated this process independently for the full phase 
dependence in our conventions. In terms of the matrix elements Zij diagonalizing the 
neutralino mass matrix ( p.lO|) : 



Z*MZ-' =Md,aa , (4.1) 
which are given functions of the parameters /i, Mi, M2, tan/3, the cross-section reads: 



0" e e 



'X.i ) 



(2 - S,,) 



<^Z + CTe + CTZe 



(4.2) 



where the three contributions to the integrated cross-section are given by 



^vrsy^ cos 



4 

w 



\Dz{s)\'{h \0 



"L\2 
ji I 



2s m,m, ^[{O-l^y]} (4.3) 



= ^^^^'^l^f 'l^f '^3^^'^-) - ' ^[(/^t)'(//r)'] ^5(me-J + (L ^ R)} (4.4) 



L\2/- fL*\2] 



-Re {miTfifO';r] him, J - s m^m, WiT fifO-l^Mm,,))} 

(4.5) 

where rrii, rrij designate the two neutralino physical masses, and the various couplings 
are defined as 

/f = V2 [(T3/ - ei) tanOw - T^i Z,^\ 
= V2Qtan^i^Z;, 

(4.6) 



U = T^i - ei sin^ 6w , Ri = -ei sin^ 6w, (4.7) 

and 

= -lz,,z*, + ^z,,z;, = o';r (4.8) 

is the coupling of the Z to neutralino i,j with O'-j^ = —O^^* (see the second ref. in [|1| 
for more details, but note a typo in eq.C77 therein.) 
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The kinematical variables are 



As 



[s — {rrii + mjY){s — (mj — m 



Dzis) 



s — rriz + imzTz 



0, 



(4.9) 



and the expressions for the integrals appearing in eqs. (|4.3|-|4-5|) are collected in the 
appendix, where we also display the relevant expressions for the amplitudes. Note that 
the complex phase factors in the mass eigenvalues (usually denoted by r/j [^-[^) are in 
our procedure automatically taken into account in eqs. ( ^4.3| -0), in both CP-conserving 
and CP-violating cases, through the redefinition 



ry 



idj /2 



(4.10) 



when the mass eigenvalue of the neutralino j species picks up a 9j phase 



Our expressions agree analytically with ref. ||3T| in the real case and with ref. |^ in 



the complex case. We obtained as well very good numerical agreement with the values of 



the cross-section displayed in ref. 23 . 



4.2 Illustrations of X1X2 production cross-section 

In fig. 1^ we plot the e+e^ X1X2 total cross-section for the same choice of parameters 
corresponding to fig. la,b, and two different choices of the selectron masses (see figure 
caption), for a LEP2 energy of 190 GeV (top figure) and a future linear collider energy of 
500 GeV. (For simplicity, we assume a negligible mixing in the selectron sector, so that 
and the physical masses). Although our input parameter choice is partic- 

ular, the behaviour of the cross-section is rather generic. As a first general observation, 
this plot illustrates the important sensitivity of the total cross-section to the variation of 
$^ and (the plots in fig. 9 are functions of but varies also for this input 
choice, according to our inversion algorithm, as is clear from the corresponding fig. lb.) 
For instance, for both choices of selectron masses illustrated in the figures, the values 
of the cross-section vary roughly by about 30-40 % when the phases are varied within 
their maximal possible range, both at LEP2 and linear collider energies. This is not too 
surprising, since although the two neutralino masses (and the two selectron masses) are 
fixed, there is an important sensitivity to the phases through the Z-neutralinos and e- 
e-neutralino couplings entering the cross-section formula, eqs. (|4.3| - |4.5|) . Now a rather 
unpleasant feature of the presence of non-zero phases is that, depending on the selectron 
masses, one may have several different $^ i^Ah) values giving the same total cross-section 
value, even though the two neutralino masses and the selectron masses are fixed (see for 
instance the dotted lines in fig. 9, corresponding to the second Mi solution at 190 GeV, 
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for mgj^,mg^ = 100, 120 GeV). These ambiguities on the phases may however be eas- 
ily resolved when looking at the corresponding values of the chargino masses (provided 
the latter are known), whose corresponding variation with $^ is illustrated in fig. la. 
This clearly stresses again the importance of looking at possible correlations between the 
neutralino and chargino sectors in general. Note also that the cross-section values in fig. 
9 are symmetric with respect to $^ ^ 27r — $^ (which also implies —^Mi, see 

e.g. fig. lb). Concerning now the Mi reconstruction ambiguity, which is more specific to 
our inversion algorithm, the plots in fig. 9 also illustrate that the total cross-section is 
generally quite sensitive to the twofold Mi solution: even though the two \Mi\ solutions 
of fig. la are very close to each others, the two solutions for the phases are not, 
thus leading to the above mentioned sensitivity. Therefore, provided the cross-section is 
sufficiently large to be measured with a reasonable accuracy [|4|, one should be able 
to resolve this ambiguity in our Mi reconstruction procedure rather easily, provided of 
course that the selectron masses are also known. 

In fig. |l^ we plot the e+e^ X1X2 total cross-section at LEP2 (top figure) and future 
linear collider energies (bottom figure), but now as a direct function of the chargino, neu- 
tralino masses and chargino mixing angle, for the same choice of parameters corresponding 
to fig. 3a,b and one choice of the selectron masses (see figure caption). The regions with 
no plots in fig 10 simply correspond to the excluded 0l range from consistency of the 
inversion, as discussed previously for the plots in figs. 3a,b. Again, the different plots il- 
lustrate the generically rather important sensitivity of the cross-section to the phases (the 
plots in fig. 10 are functions of the mixing angle 0/,, which corresponds for fixed chargino 
and neutralino masses to varying $^ and $Afi, as is clear from the corresponding fig. 3b). 
One should keep in mind that since the mixing angle 0l (more precisely cos 20/,) can be in 
principle determined from the measurement of the chargino pair producion cross-section, 
plots of the neutralino pair production like in fig. 10 may be viewed as a function of the 
chargino pair production, for fixed chargino and neutralino mass values. This illustrates 
further the strong correlations between the chargino and neutralino physical parameters. 
As for the ambiguities in our inversion algorithm, they come in this case too from the 
twofold solution in the reconstruction of Mi. Again, given the generically quite differ- 
ent cross-section values corresponding to the two Mi solutions in fig. 10, this ambiguity 
should be easily resolved through the measurement of this observable. 

5 Conclusion 

In the present paper, we have derived a purely algebraic algorithm to reconstruct the 
(unconstrained) gaugino sector Lagrangian parameters fi = |/i|e**'', M2 = IM2I and 
Ml = |Mi|e**"i, directly from the physical chargino and (some of the) neutralino masses. 
Our construction exhibits in a more direct and systematic way the non trivial correlations 
among the physical chargino and neutralino physical parameters, which exist even when 
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the maximal possible phase freedom of the unconstrained MSSM parameter space is con- 
sidered, and which may be very hidden or cumbersome to extract in the more standard 
approach of systematic scanning over the basic parameters. Our approach should be use- 
ful in particular in the case where only a subset of the minimal required input is available, 
and we illustrated with several typical such scenarios the kind of theoretical consistency 
constraints that could result, in addition to independent experimental constraints. 

We have also recalculated the neutralino pair production cross-section e+e~ XiX^ 
within our framework and phase conventions, and illustrate its sensitivity to the phases of 
/X and Ml, which should be in particular useful to resolve an ambiguity in the Mi param- 
eter reconstruction, intrinsic to our algorithm. A careful comparison of the chargino and 
neutralino pair production together with the use of the inversion algorithm should also 
give non trivial correlations among the different chargino/neutralino masses and cross- 
section observables. Those possible constraints may be studied in a more systematic way, 
taking into account also the additional direct or indirect constraints from LEP data, mea- 
surements of the electric dipole moment, etc... 



Finally, a refinement of our tree-level algebraic inversion procedure would be to include 
radiative corrections. However, these corrections are expected to be generically much 
smaller for the gaugino masses than for the chargino pair production cross-section , 
and as such, would probably affect more the reconstruction algorithm of the chargino 
parameters than that of the neutralino ones. 
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Appendix: Neutralino pair production formulas 

We collect here some formulas needed for the expression of the neutralino i,j pair pro- 
duction cross-section. The integrals in ( [4.3H4.5| ) read 

4 q^y s 



him] 



1 8(2s(-g2 + EiEj) - ^{Ei + Ej){-2m^ + mf + m|) + (-2m2 + mf + 
2} - Vs + {-2m? - ^{Ei + Ej) + mf + m^f 

2L[m](2m^ — mf — m^) 

L[m] 



2\2\ 
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L[m] 



-2m2 J^mf + m]- -/s^Ei + Ej - 2q) 



(5.1) 



The square amphtudes for the neutrahno i, j pair production involve s— , and li— 
channel exchanges. They read: 



-2s m,m, D,,{t)D,Mmfi-)Miy] + ^ ^)} (5-3) 

25i[r,(r, + r„)*] = 44{ 



Um[Dzis)f^jt:o;r]-sm.m,^Dz{s)ftjtO';t] ] + 
Z^e-,(«)[ e.,(«)3^[/^z(s)/f/fO;^]-sm,m,5R[D^(s)/f/fO;f] ] } 
-i?e{A.Jt)[ 6.(t)5J[/?z(3)/,f/ro;:^]-5m,m,$R[D;,(s)/,f/f^^^ ] + 

D.Mi iiAu)mz{s)fifh';0]^^-sm,m,mz{s)fiffi':0]^^ ] } } 

(5.4) 



and the other functions and couplings as defined in section 4. The differential cross-section 
is given by 



where 



iij{x) = (m^ -x)(m^2 -x) 



(5.5) 



(5.6) 




(|r,|2 + 2^T,{Tt + TuY] + \Tt + r„n 



(5.7) 
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Figure 1: the twofold \Mi\ (top figure) and $Mi (bottom figure) reconstruction with 
input choice Ifi] = 100 GeV, Ms = 120 GeV, M^^ = 40, Mn^ = 80; tan/3 = 2. Also 
shown are the corresponding chargino mass values. 



20 



(GeV) 



400 




(rad) 3 



2 - 



- 



-2 - 



$Mi(l) 



NO SOL. ZONE 



6 

$^ (rad) 



Figure 2: same as fig. 1, with input choice |//| = 100 GeV, M2 = 120 GeV, Mn^ 
= 60; tan/5 = 10. 
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Figure 3: M2, (top figure) and $Mi (bottom figure) full reconstruction with 

physical input choice = 80 GeV, = 180 GeV, Mn^ = 40, Mn^ = 80; tan/3 = 2. 

22 



200 

(GeV) 




(rad) 



3 


1 


1 1 M 1 1 1 - 


2 


- % 




1 













1 







2 




$M.(1) 


3 


1 


1 1 1 1 1 1 - 



0.2 0.4 0.6 0.8 1 1.2 1.4 

(pL (rad) 

Figure 4: same as fig. 3 with input choice = 80 GeV, = 180 GeV, Mn^ = 40, 
M7V2 = 80; tan/3 = 5. 
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Figure 5: Correlations between arbitrary neutralino masses from consistency of the 
inversion. Input parameters: — 80 GeV, = 180 GeV (indicated by the two 

straight hues), 0.37 < (pLirad) < 0.52; tan/3 =2. 
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Figure 6: same as fig. 5 for a difi^erent input choice: M + = 80 GeV, M + = 180 GeV, 
(f)L{rad) = 0.5 ; tan/3 = 5 
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Figure 7: Correlations in the plane from consistency of the inversion, with 

incomplete knowledge of the input parameters: — 80 GeV, 80 < < 1000, 

<4)L< 7r/2, MiVi = 40, Mn^ = 70; 1 < tan/? < 3. ' 
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Figure 8: same as fig. 7 with a different input choice: = 200 GeV, 200 < < 

1000, <(f)L< 7r/2, Mjvi = 150, Mn^ = 190; 1< tan/3 < 3.' 
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Figure 9: Xi)X2 production cross-section at ^/s — 190,500 GeV, versus for two 

different choices of selectron masses, and input choice as in fig. 1. Plain (dotted) lines 

correspond to the first (second) Mi solution of fig. 1. 
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Figure 10: Xi^ X2 production cross-section at ^/s = 190, 500 GeV, versus the chargino 
mixing angle (pL, and for the same physical input as fig. 3. Plain (dotted) lines correspond 
to the first (second) Mi solution of fig. 3. 
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